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INTRODUCTION 


In this report, three coding schemes designed for improved data communi- 
cation are evaluated. In Part A, four block codes are evaluated relative to 
a quality function, which is a function of both the amount of data rejected 
and the error rate. 

Part B is an evaluation of the Viterbi Maximum Likelihood Decoding Al- 
gorithm as a decoding procedure. This evaluation is obtained by simulating 
the system on a digital computer. 

In Part C, Short Constraint Length Rate 1/? 'Quick-Look' Codes are 
studied, and their performance is compared to general nonsystematic codes. 
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PART , A 



PERFORMANCE OF BLOCK CODES 


I . Introduction: 

Although the use of error control coding techniques in digital space com- 
munication systems has become fairly routine in recent years, there still exists 
a great deal of uncertainty as .to the actual effectiveness of coding in achieving 
more reliable communication. The reason for this is to be found in the fact that 
the commonly used performance parameters do not take Into account all the perti- 
nent aspects of the coded transmission system. Thus, for example, the widely 
used Probability of Word Error criterion totally ignores the possibility that 
the decoder may incorporate some degree of data rejection. Likewise, the mini- 
mum distance criterion, another popular measure of code performance, is com- 
pletely independent of the decoding algorithm and several other important system 
factors. 

As a consequence of this state of affairs, it is virtually impossible to 
compare, say, a coding system with error correction and data rejection to one 
with error correction alone, using any of the existing criteria of performance, 
and It Is therefore of value to define and evaluate measures which incorporate 
most, if not all, of the quantities affecting the overall system reliability. 

This is the objective of the present work. 

II. Definition of Performance Measure: 

For the simple types of block codes normally employed in space communi- 
cation systems, the complexity of the encoder and decoder is of little conse- 
quence, since the use of integrated circuit technology allows the construction 
of the basic components In an inexpensive fashion. Furthermore, the complexity 
is essentially independent of the particular code-decoder used. 
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The processing speed is generally a function of the type of logic used and 
the technology in the construction of the integrated circuits. Although one 
could probably obtain cost figures as a function of processing speed, the im- 
portance of these costs in the overall system considerations is difficult to 
assess. Also, as with complexity, processing speed is not a strong function 
of the code-decoder combination. 

Thus, the important factors determining the overall coding system perform- 
ance are: 

1. The accuracy of the data after decoding, 

2. The amount of data rejected by the decoder, 

3. The amount of redundancy in the code, and 

4. The relative importance of data accuracy, data 
rejection, and data transmission rate. 

Let us consider a situation in which N blocks of received digits from a 
binary (n , k) block code are to be decoded. The decoder generally rejects N-X 
blocks, leaving X blocks after decoding, of which V are correct. (See Figure 1) 


4 N-X BLOCKS REJECTED 


N BLOCKS 


DECODER 


X BLOCKS 

Y CORRECT BLOCKS 


FIGURE 1. GENERAL DECODER CONFIGURATION 


The amount of data passed by the decoder is measured by the quantity 


F, = n E{X} ’ 
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the accuracy of the data after decoding is measured by the quantity 
F _ E{Y} 

2 “ eTxT ’ 

and the amount of redundancy in the code is measured by the quantity 

p J< _ number of data digits per block 

3 n total number of digits per block * 

Here E is the usual expectation operator. 

We also define a quantity 0 < a < 1 which measures the relative importance 
of data accuracy and data rejection. 

As an overall measure of performance of the code-decoder combination, we 
then take quanti ty 

F = 1 - 

as a function of the energy per information bit-to-noise ratto, E^/N o . 

When the N blocks are transmitted independently of each other and are 
treated as such by the decoder, 1 - F reduces to the probability of word rejec- 
tion. For a decoder with no data rejection, F^ becomes the probability of cor- 
rect decoding. Thus, in the two limiting cases a = 0 and a = 1, F reduces to 
the probability of word error and word rejection, respectively. 

III. Evaluation of F for Hamming Type Block Codes Over the Binary Symmetric 
Channel : 

We assume that N blocks are transmitted independently and with equal pro - 
bability over a binary symmetric channel whose digit error probability is 
p = 1 - q. The codes of interest are of two types: The standard (n, k) Hamming 

code described by the parity check matrix 
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whose columns are all 2 - 1 nonzero binary m-tuples (m any integer greater 

than 2), and a modified Hamming code whose parity check matrix differs from 
the above only in having an additional row of ones on top. The first code has 
block length n - 2 m - 1, k = n-m information digits and minimum distance 3 
and is thus able to correct all single errors. The second code has the same 
block length, k ” n ~ m - 1 and minimum distance h and can be decoded in either 
a single-error-correcting, double-error-detecting mode or a t r i pie-error- 
detecting mode* 

For both codes, the first step in decoding a received block v = (v^ , v^, 

• «■ • > v n ) consists of determining its syndrome* This is a binary (n-k)-tuple 
given by 

s = vH T 

where T denotes matrix transposition and the multiplication and addition opera- 
tions are modulo 2. 

We now consider four cases, including, for purposes of comparison, the 
uncoded transmission of data blocks of length n. 

Case 1. No Coding - (n, n) Code 

Decoding Rule: Pass every block unchanged 

Evidently, X = E{X} = N and since a block is correct at the 
decoder output if and only if it is correct at the input, we 
have E{Y} = N n 

q 
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Hence = 1 ; = q n 

and 


F - 1 - q n(, -“> 

Since F^ = 1 , the relation between the channel 

p and £, /N is p = i - erf E, /N 
r b o 2 bo 


error p robabl 1 1 ty 


Case 2. Single-Error-Correcting Hamming Code 

Decoding Rule: If the syndrome is zero, pass the block, tf 

the syndrome is not equal to zero, assume a single error has 
occurred, correct it, and then pass the block. 

Again, E{X} = N. For E{ Y } we have 

E{Y} = N{ Probabi 1 i ty that a block has no error or a 
single error before decoding} 

= N{q r ' + nq n - ^ p} 


Therefore, F^ = — , F 2 = ^ + nd" 1 P 

and 


F - 


i / n . n-1 

1 “ (q + npq 


) 


where p 



E b /N 


o 


Case 3. Single-Error-Correcting, Double-Error-Detecting Hamming Code 
Decoding Rule: If the syndrome is zero, pass the block. If 

the first digit and at least one of the remaining digits in the 
syndrome are one, assume a single error has occurred, correct 
it, and then pass the block. For all other syndromes, reject 
the block. 


We have F, = , 

3 n ’ 
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E{X} 

Fj = — ^ — = {Probability that a block has zero syndrome 




or the first and at least one of the remaining 
digits equal one} 


n-2 i 2 i 


+ [ (2 i n + 1) - A 


n-2 1 +1 2i + 1 


2i + 1- 


where Aj is the number of codewords of weight j of the Single- 
Error-Correcting Hamming Code, 
and for we obtain 


E{ Y} 1 

^2 = eTxT = F~” ^ ro ^ a bi I i ty that a received block is correct 
or has a single error} 


Case 4. Triple-Error-Detecting Hamming Code 

Decoding Rule: If the syndrome is zero, pass the block. Other- 

wise, reject the block. 

u c k- 1 

Here, F 3 - — , 

E{ X} 

F.j = — — = {Probability that a block has zero syndrome} 


n-1 


2 



n - 2 i 

<1 


2 i 
P 


> 


and 


F 


_ E{Y} _ £ 

2 “ eTxT " F 2 


The Hamming code weight spectra required for Cases 3 and 4 may be obtained 
as the coefficients of the polynomial. 
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f <*> ■ nTf 


n-1 n+1 

{ (l+x) n + n(l+x) 2 (1-x) 2 ) 

where A. is the coefficient of x 1 . 

IV. Evaluating the Performance Criterion: 

A Fortran language program, reproduced in Appendix A, was written to evaluate 
the function F for the four cases described above. The program calculates F for 
101 equally spaced values of E b /N Q ranging from 2 db to 10 db and all values of 
redundancy from m = 3 to m = 10. 

A major part of the program is devoted to calculating the coefficients of 
the function 

n~ 1 n+1 

f (X) « ((1+X) n + n (1+X) 2 (1-X) 2 ) 

which are used in Cases 3 and 4. The main difficulty in this computation is 
that some of them have magnitudes on the order of 10^^ for large values of n. 
Overflow on the IBM 370 occurs with numbers as small as 10 . To overcome this 
problem, most calculations are done using logarithms. Thus, for example, the 
logarithms of the coefficients of (l+X) n are stored in an array called LGC0EF. 
Similarly, the components of (1+X) n ^ 2 and (1+X) n+ ^ 2 are stored in LCNM1 and 
LCNP1 , respectively. Note that LCNP1 contains the logarithms of the coeffi- 
cients of (1+X) n+1/2 and not (1-X) n+1/2 , since the latter has negative coeffi- 
cients whose logarithms do not exist. The variable SIGN, which always equals 
±1 , is used to convert the coefficients of (1+X) n to the coefficients of 
(l-X) n+ ^^ when making calculations of f (X) . 

A special procedure is used throughout the program to achieve addition 
of these very large numbers. Obviously this addition cannot be achieved di- 
rectly using logarithms. To illustrate this procedure consider adding the 
numbers A = 7-3147 x 10 298 and B » 2.1532 x 10 295 given the logarithms of these 


numbers. 
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ALOG = log (A) = 298.864196 « 3.864196 + 295 . 

BLOG = 1og(B) = 295.333084 = 0.333084 + 295- 
Let Z = X+Y and ZLOG = log (Z) . 

ZLOG = logOO 3 - 8 ^ 196 + 10 °-233084 ) + 295 

= log (7314. 7 + 2 . 1532 ) + 295 
= 1og(73l6.8532). + 295 
= 3.864324 + 295 
= 298.864324 

Thus, the log of the sum has been calculated using numbers no bigger than 
7316.8532. Since the calculations on the IBM 370 have only 16 significant 
figures, numbers whose magnitudes differ by more than lO 1 ^ are not added by 
the above method. In this case the sum is set equal to the larger of the 
two numbers. 

The coefficients of f(x) are calculated using the aforementioned tech- 
niques and stored in an array called RIGHT. The variable RINOM is set equal 
to the logs of certain binomial coefficients, and it is used in calculating 
terms of FI CAS 3 (Fj for Case 3; i.e., F ] for SEC-DED) of the form 

«f> - Aj^'V. 

These terms are stored in C0EFC3. 

In order to calculate the parameters for each case for any particular 
value of S/N, values for q and p, which are dependent on the code rate, must 
be evaluated. The dependence on code rate requires calculations of Q1 and PI, 

Q 2 and P2, and Q3 and P3 for use with Cases 1, 2 , and 3 and 4 respectively. 

Since Case 4 has the same code rate as Case 3, Q3 and P3 are applicable to both. 

IV. Numerical Results and Conclusions: 

In Figures 2-49, we have plotted the performance measure F as a function 
of the signal-to-nolse ratio E^/N o of the binary symmetric channel in db, for 
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all four cases described above, and for n = 7, 15, 31, 63, 127 , 255, 511 , 1023 , 
and a — 0, 0.1, 0.3, 0.5, 0.8, 1.0. We draw the following conclusions. For 
the extreme values of alpha the relative importance of each case remains fixed 
for all values of E^/N^ When a = 0 , indicating an interest only in the quality 
of the output, the relative ratings from best to worst are: TED, SEC-DED, SEC, 
No coding. As expected, when a - 1, indicating an interest only in the quantity 
of output data, the relative ratings are just opposite to the a = 0 cases. For 
a = .5 and n £ 15, the relative order also remains fixed: SEC, SEC-DED, No 

coding, TED. Notice that the extreme cases of large quantity of output achiev- 
able with no coding and high quality of output achievable with a TEC system 
are both given poor relative ratings for this value of alpha showing no pre- 
ference of quantity over quality or vice-versa. Also, for a = .5 n = 7 , no 
coding becomes preferable to SEC-DED at s i gna 1 - to-noi se ratios below 4.1 db. 

This would be due to the increased data rejection by a SEC-DED decoder as the 
noise becomes greater. 

A preference for quality over quantity without total disinterest in the 
latter is explored by the a = .1 and a = .3 cases. No coding for these values 
of alpha is never preferable to any other of the choices considered except TED. 

As S/N increases TED becomes relatively less desirable as triple errors become 
less likely and its low transmission rate becomes dominant. Similarly, as 
E b /N o increases SEC-DED becomes less important than just SEC. However, as the 
block size increases, TED and SEC-DED become more important since the proba- 
bilities of the errors these decoders are designed to correct increase. 

Finally, when quantity is somewhat preferred over quality as with a = . 8 , 

as might be expected the relatively extreme quality achieved by TED is shown 

to be undesirable for all values of S/N tested since this quality is achieved 

at the expense of quantity. The SEC decoder which employs no data rejection 

yet achieves some degree of error correction is found to be the best of all 
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four cases for all values of E,/N considered. At very low noise levels the 

DO 

error correcting properties of SEC-DED make it more desirable than no coding, 
while at high noise levels the data-reject ion of SEC-DED make it less de- 
sirable than no coding. For example, for n = 1023, SEC-DED is preferable to 
no coding for E b /N Q from 4.4 to 6.9 db. 

More insight into the nature of this function can be gained by looking 
at what is necessary to achieve a desired level of effectiveness. A typical 
example is shown in Table 1. Here the desired value of F is set at .01. 

With no coding or SEC, a higher value of E b /N^ is required as quality becomes 

more preferable. However, with SEC-DED or TED a lower value of E./N is re- 

b o 

quired to achieve the same value of F as emphasis is switched to quantity. 


E./N !N DB REQUIRED TO 

D O 

ACHIEVE F = .01 FOR n = 255 


Alpha 

No Coding 

SEC 

SEC-DED 

TED 

.1 

8.88 

7.36 

6.64 

7.68 

.3 

8.72 

7.20 

6.88 

8.40 

• 5 

8.56 

7.04 

7.12 

8.72 

.8 

8.00 

6.64 

7.28 

8.96 


TABLE 1 
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DIGITAL SIMULATION OF THE VITERBI 
MAXIMUM LIKELIHOOD DECODING ALGORITHM 

I . I nt roduct i on : 

The Viterbi algorithm is a method for determining the most likely sequence 
of states of a time-discreet Markov process; and, as such it is an optimal 
method for decoding convolutional codes. An evaluation of the effectiveness 
of this algorithm as a decoding method is accomplished herein through simula- 
tion on an IBM 370 computer using a main program written in the Fortran lan- 
guage and three subroutines written in Assembler language. 

II. The Simulation Procedure: 

A block diagram of the simulation is shown in Figure 50 . A pseudo-random 
number generator is used to independently generate binary source bits of equal 
probability and binary noise bits whose probabilities depend on the assumed 
channel characteristics. The source bits are encoded with the appropriate 
parity check bits in blocks of two. Each bit is added to a noise bit using 
modulo-2 arithmetic (simulating channel noise). The information noise bits 
and parity check noise bits are generated independently. The corrupted bits 
are then decoded using the Viterbi algorithm. Accuracy of the decoding al- 
gorithm is measured as 

E - limit E = (Number of correct information bits)/(Total bits) 

n ■+■ co 

where E n is the ratio after information bits have 
been transmi tted. 

III. The Viterbi Algorithm: 

Given an observed output sequence Z = , z 2 , . . . , * ) , the purpose of 

the Viterbi algorithm is to determine the most likely input sequence X = (x Q , 
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BLOCK DIAGRAM OF 


THE VITERBI ALGORITHM SIMULATION 


FIGURE 50 






» • • • , x k ) • The subscripts refer to the discreet time states. Since the 
process is assumed to be Markov, the probability of state x. + .j depends only on 
the state x. : 

i 

i.e., Pr (x. + j j X Q , Xj , . . . , x.) — Pr (x | x j ) . 

The channel is assumed to be memoryiess so that the observed output z. at time 
i depends only on the transition from state x. to state x. + 1 * This transition 
is symbolized as t.. We want to determine the maximum a posteriori Pr(X, Z) . 
Because of the aforementioned Markov and memoryless assumptions: 

Pr(X, Z) = Pr(x)Pr(z|x) 

k-1 k-1 

- " P r ( x j,i » x.) 7T Pr (z. 1 t.) 
i=0 1+1 ' i=0 11 

The Viterbi algorithm is a method of determining the shortest path be- 
tween two points. We, therefore, assign a 'path length* between each pair 
of possible states from time = i to time = i + 1. This length lambda(t.) is 
defined as 

lambda (tj ) = -In Pr(x f+1 |x } ) - In PrUjt.). 

The total length for some input sequence X would be 

k-1 

-In Pr (X, Z) = E lambda(t.) 

1=0 1 

Since path length is a negative logarithm of the probability, the shortest 
(critical) path length between two points (i.e., the initial and final states) 
would be the one with the highest probability. This is the maximum a posteriori 
probability we are seeking. The Viterbi method of finding this critical path 
Is based on the observation that at any given time i, each state x. has associated 
with it a shortest path to the Initial state. This shortest sequence is called 

A 

a survivor, designated X(x.). The path length of survivor X(x.) is designated 


62 



gamma (x.). Extending these survivors to time i + 1 requires merely adding 
the appropriate digit ("bit" for our purposes) to the existing survivor and 
adding the corresponding path length gamma (x. +1 , x.) for comparison purposes 
in determining the survivors for each state at time 1+1. At the end of the 
sequence (time = k) the survivor corresponding to the state with the shortest 
survivor path length is optimal. 

For the purpose of convolutional code decoding* the states correspond to 
the possible binary state permutations of a block of m shift registers. As- 
suming that for the source Pr (0) = Pr(l) = i, it follows that for all possible 
transitions between states the term Pr(x [+t |x.) = i; and, since it is a con- 
stant for all possible transitions, it may be ignored when calculating the 
optimal path. Thus, only the term Pr(z.|t.) is significant. For systematic 
codes, the observation z. corresponds to both the information bit and the parity 
bit received as a block at time = i. For non-systematic codes, z. corresponds 
to a block containing a parity bit for each subgenerator polynomial. These 
probabilities are pre-calculated for each state and each possible received 
block before decoding begins. 

Since we are concerned with a communications system with a semi -inf ini te 
number of bits transmitted, corresponding to a semi-infinite sequence, and 
since storing the resulting semi -inf ini te survivors is impractical, a limit 
must be placed on the number of bits stored as a survivor (i.e., the survivor 
length). Call this limit delta. Thus, at time = i, a decision must be made 
concerning the bit at time = i - delta (i minus delta). This survivor trans- 
action becomes insignificant for delta large enough because survivors tend 
to converge to the same state nodes. 
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!V. The Simulation Program for Rate i Codes: 

A Fortran language main program was used In conjunction with three custom 
written Assembler language subroutines. The main thing to be aware of when 
using the program is that the delta defined in the program is one greater than 
the corresponding delta as defined in the Viterbi algorithm {e.g., if you wish 
to get results for delta = 75, set delta = 76 in this program!). For each 
state at time i there are two possible states to which it can branch at time 
i + 1 (one of which has an incoming 0 bit, the other has an incoming 1 bit). 
These possible transitions are stored in an array called NEXT. NEXT (1,1) cor- 
responds to the branch of state I with an incoming 0; whereas, NEXT(I,2) cor- 
responds to an incoming 1. 

Probabilities which determine the path lengths are calculated prior to 
the main iterative loop. These calculations are done for all four possible 
two bit permutations corresponding to a received block containing an informa- 
tion bit and a parity check bit in the systematic case, or two parity check 
bits in the non-systemat i c case. Array POFZLN stores these predetermined 
path lengths. Thus in the main iterative loop the increase in the total path 
length GAMMA of each state can be determined by a simple table reference 
(i.e., POFZLN). Survivors are stored and saved by arrays SURVIV and SAVUR, 
while the corresponding path lengths are saved using arrays GAMMA and SAVE. 

The appropriate subgenerator vectors are stored in the array GEN. In 
the case of a systematic code the second subgenerator is a 1 followed by m - 1 
zeros. For example, the subgenerators for an m = 5 systematic code would be: 
10000 and 11011. Note that the subgenerator 10000 merely generates the in- 
formation bit. SHIFT2 saves the contents of the simulated encoding shift 
reg i sters . 
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All random bit generation is done independently for each application. 

Source bits of equal probability are generated using a Scientific Subroutine 
Package member called RANDU. These source bits are stored in PRSOUR. The 
generated noise bits (Pr(0) - q) are stored in PNOISE. Decoded bits are stored 
in PROUT for comparison with the original source bits in PRSOUR. 

Many different counters are used to keep track of time states corresponding 
to source bits, noise bits, and output bits, I0UT determines the printed in - 
crements for E^, During the course of this research, I0UT was set equal to 
1000 so that the accuracy was printed out for n = 1000, 2000, 3000, etc. 

The n refers to the number of decoded bits and is called NDEC0D within the 
program. 

Read and punch statements are included to save the Information necessary 
to restart the program where it left off. This feature is desirable to enable 
the programmer to check the convergence of the accuracy figures and compare 
them with other data in order to determine the desirability of decoding a 
greater number of bits. However, due to the fact that hexadecimal double 
precision accuracy used by the IBM 370 computer is equivalent to about 16.7 
decimal digits, and the data cards are punched with decimal numbers, there is 
a slight loss in accuracy that is sometimes noticeable but generally insignificant. 

Three Assembler languages subroutines were written to expedite the execu- 
tion of the program. These are C0PYAR, SHIFT, and TESTBT. C0PYAR simply copies 
SAVSUR into SURVIV. SHIFT is used to shift the survivor of a row in SURVIV, 
add a 0 or a 1, and transfer the resulting survivor to the row specified by 
NEXT in the array SAVSUR for the next time increment. TESTBT determines whether 
the bit at time = k delta (Viterbi definition of delta) is a 0 or a 1 in the 
survivor of the current state whose total path length is the least. 
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Typical decoding rates are shown in Table 2. The number of bits refers to 
the total number of both parity and information bits. To obtain the number of 
information bits decoded per second, multiply the rates shown in Table 2 by the 
code rate, which in all cases explored here is i. Note that decoding an m = 7 
code is approximately twice as slow as decoding an m = 6 code since the latter 
has half the possible states of the former. 

TYPICAL DECODING RATES 


m 

del ta 

b i ts/ sec 

5 

31 

578 

5 

52 

506 

6 

31 

32k 

6 

59 

275 

7 

27 

165 

7 

59 

145 


TABLE 2 


V. Program for Simulating A Viterbi Decoder 
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'NXMTD = 0?^ 


+ ^ Yes 


Initialize: ERRORS, 

NTOTAL, NCOUNT , 
NDECOD , PARCK(N), 
SHIFT2 (M) , 

GAMMA (NSTATE) , 

URVIV (NSTATE, NCOL)_. 


Read : ERRORS , 

NTOTAL, 

PNOISE (NCOUNT, N) , 
PRSOUR (NCOUNT) , 
PARCK (N) , 

GAMMA (NSTATE) 

SURVIV (IDEPCP ,NCOL) , 
SHIFT2 (M) . 


Define IOUT. 
Initialize NPLOT. 
Print STON and Q. 
Establish P with a 
lower limit of 
l.D-75. 

Define: QLN, PLN. 


Calculate POFZLN. 

This array contains the 
predetermined branch 
lengths based on the 
probability that a 
given state is correct 
for the block of bits 
received, i.e. , 

POFZLN (I, LOC) contains 
the probability, ex- 
pressed as a path length 
that state I is correct 
for block number LOC 
which can be any of 2**N 
ermutations. 
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Call COPYAR to copy 
SAVSUR into SURVIV. 


NDECOD< 0? 


Call TESTBT to make 
a final decision on 
the bit at time 
NXMTD- (DELTA-1) 
based on survivor, 
of state IMIN. 

Store this bit 
in PRQUT. 


NDECOD 

>IOUT 


Increment NDECOD and 


IXMTD>MAXDEC? 


Set IOUT = 
NDECOD- 1 


Increment NPLOT 


Determine cumulative 
number of noise bits 
on each block bit 
corresponding to 
NTOTAL information 
bits using PNOISE 
nd store in PARCK, 


me total 
e decoded 
parincr PR 
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PROGRAM FOR SIMULATING A VITERBI DECODER 


FOR STORAGE REASONS DELTA IN THIS PROGRAM SHOULD BE SET TO ONE MORE 
THAN THE CASE BEING SIMULATED. 

I.E. TO SIMULATE DELTA=I3, SET DELTA=14. 

THE FOLLOWING MINIMUM ARRAY DIMENSIONS MUST BE ALLOCATED: 

PNOISE I lOUT+OELTA,N) 

PRSOUR ( IOUT+DELTA I 
PROUT (I OUT) 

PNPLOT ( I OUT , N ) 

GEN (M,N) 

SHIFT2 (M) 

L (M) 

NEXT INSTATE, 2) 

PAR IN) 

OUTPUT INI 

XMTPLT INUMBER OF LINES PRINTED 2 1 NUMBER OF INFORMATION 

BITS DECODED THIS RUN/IOUT)) 

FPLOT (SAME AS XMTPLT) 

PARCK IN) 

POFZLN INSTATE »N2=2**N ) 

GAMMA INSTATE) 

SAVE INSTATE) 

IXNOIS IN) 


SURVIV AND SAVSUR MUST BE DIMENSIONED EXACTLY (NROWS.NCOL) WHERE 
NROWS MUST BE AT LEAST 2**M AND NCOL MUST BE AT LEAST IDEPTH 


INTEGER PARSUM 

INTEGER*2 PNO I SE 1 1064, 3 ) » PRSOIJR 1 1064 ) 

INTEGER GAMCHK, TEST 
INTEGERS PROUTUGOD) 

INTEGER PNPLDTI IDOGt 3) 

INTEGER*2 GEN I 8 , 3 ) T SHI FT2 I 8 > , JKL , IOUT , M, K2 , OELT A , DELTA1 

INTEGER*2 DELMi , NSTATE ,L I 8 ) , LR , NEXT I 2 56 , 2 ) , P ARt 3 ) .PARITY 
INTEGER*2 OUTPUT I 3 ), LOC ♦ SOURCE 
INTEGER SURVIVI 2 56 , 2 ) , SAVSUR I 256 , 2 ) 

INTEGER XMTPLT! 1000) 

INTEGER PARCK I 3) 

REAI -*8 POFZLN I 256,4 I , GAMMA (256) , SAVE! 256) 

REAL*8 FPLOT ( 100 ) 

INTEGER NTOTAL, ERRORS , NPLOT, NCOUNT, NDECOD, NXMTD , MAXB I T, IXNOIS (3) 
REALMS Q»P, DFLOAT, $TON»OERF, DSQRT»QLN, PLN 
REALMS DLOG,GAMMIN»GAMMAX 

C MNPLOT AND PNPLOT STORE INFO AND PARITY NOISE COUNTS RESPECTIVELY 

C MAXB I T IS THE MAXIMUM NUMBER OF BITS TO BE XMTD 

READ1104,MAXB IT 
1 104 FORMAT! 17) 

C N IS BLOCK LENGTH 

REAOllOD.N 
N2=2**N 

C K IS THE NUMBER OF INFORMATION BITS REPRESENTED BY A BLOCK 
C OF LENGTH N. 

C THIS PROGRAM WILL ONLY SIMULATE CODES OF RATE 1/N. 

C THEREFORE K MUST ALWAYS EQUAL 1, 
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K= 1 

K2=2**K 

C M IS CONSTRAINT LENGTH 
READ1 100,M 
1100 FORMAT! ID 

READ11C4,NR0WS 
READ1 If, A . NCOL 

C GEN IS SURGENERATOR VECTOR 
DO 1506 J= It N 

1506 READ1102, (GEN(IfJ)f I=1,M) 

1102 FORMAT (911) 

PRINT 608 

698 FORMAT! • I'D •> 

C OELTA=SURV IVOR LENGTH PLUS ONE 
READ1106, DELTA 
1106 FORMAT <13 1 

C ST ON IS THE S IGNAL-TO-NOISE RATIO IN DECIBELS. 

READ11D8,ST0N 
1108 FORMAT! D15. 7) 

Q=.5O0+.5OC*OERFlDSQRT(DFLOAT{K> /OFLOAT(N)*10.DO**($TON/10.DC ) ) ) 
C . NXMTO=TOTAL SOURCE BITS XMTD 
RE AD?'“-i n * NXMTD 

2000 FORMAT! 120) 

C IXSOUR AND IXNQIS ARE RANOOM INTEGERS USED AS STARTING 
C POINTS FOR THE RANDOM NUMBER GENERATOR. 

READ2 000, IXSOUR 
DO 1510 1=1, N 

1510 PEA02005' , IXNOIS ( I ) 

UELM1=DELTA-1 

C IOEPTH IS THE NUMBER OF WORDS REQUIRED TO STORE ONE SURVIVOR. 
I0EPTH=DELMl/32+l 
IDEPCP=NC0L+1-I DEPTH 
I8IT=DELTA-32* ( IDEPTH-1) 

MAXDEC=MAXBIT+DELM1 

NSTATE=2**M 

M1=M+1 

DO 9 I = 1 , M 
9 L ( I ) = 1 

DO 16 1=1, NS TATE 

00 II LR=1,M 

1 F ( L ( LR I ) 1 1 , 14, 11 

11 H LR ) = 0 

GO TO 15 

14 L(L«)=1 

C NEXT IS AN ARRAY CONTAINING THE NUMRERS CORRESPONDING 
C TO THE 2 STATES THAT ANY GIVEN STATE CAN BRANCH TO. 

15 NEXT ( I, 1 )= 1 
DO 30 LM=2,M 

30 NEXT! 1 ,1)=NEXT( I , 1 > +2**( LM-2 ) *L < LM) 

16 NEXT! I,2)=NEXT(I,1>+2**(M-1> 

PRINT1138, IXSOUR 

1138 FORMAT! IXSOUR=* , 120) 

DO 1139 1=1, N 

1139 PRINT1140, I , IXNOIS! I ) 

1140 FORMAT! *0' , ' IXNOIS! • , II t ' >* » » 120 ) 

IF(NXMTD.EQ.O)GO TO 1300 
RE AO 2000, ERRORS 
R6AD2000 , NTOT AL 
NC0UNT=0ELM1 
NDECOO= 1 

2001 FORMAT (7211) 
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DO 1550 J=1 iN 

1550 READ2001 f ( PNOISEI I , J ) « 1 = 1 v NCOUNT > 

READ2001 9 \ PRSOUR (1) ? ]*lf NCOUNT ) 

DO 1551 I=liN 

1551 READZOOOt PARCK ( 1) 

DO 13 )2 I = lf NST ATE 

1302 REA D2 002 , GAMMA I 1 ) 

20C2 FORMAT ( D23 * 16 ) 

DO 1303 1=1 tNSTATE 

DQ 13C3 J = I DEPCP j NCOL 

1303 READ 2000 f SURVIVI I , J ) 

READ200 l 9 ( SHIFT2( I) , 

GO TO 1301 

1300 CONTINUE 

C INITIALIZE RANDU SOURCE AND NOISE 
C ERRORS-NUMBER OF DECODING ERRORS 
ERROR$=0 

C NTOTAL=TOTAL BITS PRINTED 
NTOT AL= 1-DELTA 

C NCOUNT PLACES OUTPUT BITS IN CORRECT VECTOR POSITION 
NCOUNT =0 

C NDECOQ=NUMB E R OF BITS DECODED 
NDECQD=2-DELTA 

C PARCK = TOT AL NUMBER OF PRINTED PARITY CHECK BITS CORRUPTED BY NOISE 
00 1552 I- 1 » N 

1552 P ARCK ( I ) =0 

C SHIFT2=C0NTENTS OF ENCODING SHIFT REGISTER 
DO 1 I = L*M 
i SH I FT 2 ( I ) =0 

GAMMA! 11=0*00 
DO 60 1 = 2, NSTATE 

60 GAMMA! !!*1.D40 

DO 61 1 = 1, NROWS 

DO 61 J = 1 1 NCOL 

61 SURV I V ! I » J ) =D 

1301 CONTINUE 

C IOUT=NUMBER OF BITS PRINTED PER LINE 
IOUT=1QDO 

C NPLOT COUNTS NUMBER OF TIMES PRINTING ALGORITHM IS USED 
NPLOT =0 
PR I NT998 * STON 

99B FORMATl , »S/N =*#D15*7> 

PR I NT704 1 Q 

704 FORMAT I t + , f T30f f 0= t f D15#7) 

P= 1 • DQ-Q 

IFCP*LT*1* D-75 ) P= 1 • 0^75 
GLN=-1.DO*DLOG<0> 

PLN=-1 • DO*OLGG t P ) 

DQ AO 1*1, M 
49 LI I ) = 1 

00 50 1*1, NSTATE 

DO 51 LR=1 t M 
I F I L ( LR M 51 • 52 9 51 

51 L ( LR ) =0 
GO TO 53 

52 L I LR ) = l 

53 DO 1500 IBITNR»ltN 

PARSUM =0 

DO 20 11*1 1 M 

I M1=M1- I 1 

20 PARSUM=PARSUM*-GEN(Tlt IBITNR >*L< IM1) 
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1500 PARI I BITNR )=PARSUM-PAR$UM fZ*2 
00 21 13=1, N 

21 OUTPUT ( 13 ) = 1 

DO 50 J=1,N2 

00 55 14=1 *N 

14MIN1=N+1-I4 

I F ( OUTPUT I I4MI N1 ) >55,56,55 

55 OUTPUT! I4M1N1)=0 
GO TO 57 

56 OUTPUT ( 14MIN1)=1 

57 POFZLN(I,J)=O.DO 

00 50 IL77=1,N 

IF ( OUTPUT ( I L77 ) .EQ.PARI IL77) JPQFZLNI I,J )=POFZLN< I , J)+QLN 

50 IF! OUTPUT ( IL77) .NE.PARI IL77DP0FZLNC I , J >=POFZLN ( I , J I+PLN 

699 CALL RANOUI I XSOUR , I YSOUR, YSOUR ) 

IX$OUR=IYS0UR 

S0URC£=0 

IF(YSQUR-0. 5)1511,1512, 1512 
1512 SOURCE= 1 

1511 NXMTD=NXMTD+1 

NCOUNT =NC0UNT+1 
PRSOUR ( NCOUNT )=SOURCE 
DO 2 1=2, M 

ll=M+2-I 

2 SHI FT2 I 1 1 ) = SHI FT 2 ( 1 1-1) 

SHIFT2I l)=SOURCE 
C PAR 1 T Y=P ARI TV CHECK DIGIT 
DO 1516 I B I TNR=1 , N 
PARSUM=0 
DO 1517 11=1, M 

1517 PARSUM=PARSUM+GEN< II, IB ITNR ) *SH I FT2 I ID 
1516 PARI IBITNR)=PARSUM-PARSUM/2*2 
DO 1530 1=1 »N 

PNOISEI NCOUNT, 11=0 
CALL RANOUI IXNOIS(I), IYNOIS,YNOIS ) 

IXNOISII)=IYNOIS 
IFIQ-YNOIS) 1531,1531, 1530 
1531 PNOISEINCOUNT, D=1 
KPAR=PAR( I ) 

IF(KPAR.EQ.O)PARI I ) = 1 
IFIKPAR.EQ.l )PAR( I 1 = 0 
1530 CONTINUE 

C LOC= LOCATION IN PROBABILITY MATRIX CORRESPONDING To BLOCK RECEIVED 
L0C=1 

00 1504 L0CSUM=1»N 
504 L0C=L0C+2** I N-LOCSUM ) *PAR (LOCSUM) 


STATES ARE NUMBERED SUCH THAT 1 AND 2, 3 AND 4, ETC. ARE 
PAIRS THAT BRANCH TO THE SAME STATES. 

I.E. OX AND IX BOTH BRANCH TO XD AND XI, WHERE X REPRESENTS 
A PERMUTATION OF M-l BITS. 


GAMCHK IS THE STATE OF GIVEN PAIR OF STATES WHICH HAS THE SHORTEST 
TOTAL PATH LENGTH GAMMA. 


DO 100 1=1 , NSTATE, 2 



GAMCHK-I 

I F ( GAMMA ( I ) .GT .GAMMA ( I + 1 > )GAMCHK=I+1 
DC 101 J=1 »K2 

101 SAVE ( NEXT ( I » J ) )=GAMMA(GAMCHK)+P0FZLN(NEXT! t , J) ,L0C) 

irowc=nexti 1,11 

IR0WI=NEXT< 1,2) 

CALL SHIFT(SURVlVfSAV$UR,NROWS* NCOL *GAMCHK,J R0WC- , I ROW 1 ) 
100 CONTINUE 

DO 139 1 = 1 ,NST ATE 

139 GAMMAU)=SAV=(I) 

L42 GAMM I N = 1 . D70 

GAMM AX=-1 • 070 
I M I N = 1 

DO 140 I = 1 , NSTATE 

1 F ( GAMM A ( I ) .GT.GAMMAX)GAMMAX=GAMMAl I ) 

If (GAMMA! D.GE.GAMMINJGO TO 140 
GAMMIN-GAMMAI I ) 

I MIN* I 

140 CONTINUE 

155 IFIGAMMAX.LE. 1 .D50)C0 TO 150 

00 141 1=1, NSTATE 

141 GAMMA! I )=GAMMA! I ) -1 . D5 1 " 1 - 
GO TO 142 

150 CALL COPYAR<SAV$UR,SURVIV,NROW$,i\lCOL) 

IFtNDECOU.LE.OJGO TO 715 
CALL TESTBTI SURVIVI I M IN , IDEPCP > , IB ! T, TE ST ) 
PROUT(NOECOD)=TESr 
715 NDECOn=NDECOD+l 

NT0TAL=NT0TAL+1 
IFINXMTO.GE.MAXDEOGO TO 750 
I F I NOECUD. GT . TOUT ) GO TO 700 
GO TO 699 
730 CONTINUE 

NPLOT =NPL0T*1 
00 776 J= 1 , N 
00 776 1=1, IOUT 

776 IF! PNOI SE! I » J I »EQ ♦ 1 ) P ARCK ! J ) =P ARCK ( J I + 1 
00 7.70 1 = 1, IOUT 

770 IF ( PR SOUR ( I ) .NE.PROUTt l ) ) £RR0RS=ERR0RS+1 

FPLOT!NPLDT)=OFLOAT!NTOTAL-ERRORS) /DFLOAT! NTOTAL ) 
DO 1561 1=1 , N 

1561 PNPLOTINPLOT, I )=PARCK ( I ) 

XMTPLT (NPLOT)=NTOTAL 
IFINXMTO.GE.MAXDEOGO TO 900 
DO 778 J= 1 , N 
DO 778 I =1 ,DELM1 

778 PMOISEI I* J )=PNQISE( IOUT+I, J) 

DO 703 1=1 , DELMi 

703 PR50UR I ! ) = PR SOUR C I OUT + 1 ) 

NCOUNT = DELM1 
ND£COD=l 
GO TO 699 

750 10UT=NDEC0D-1 

GO TO 700 
900 CONTINUE 

PR I NT 500 , M 

500 FORMAT! •-*, 'CONSTRAINT L£NGTH=«,I3> 

DO 1565 J= 1 , N 

1565 PRINT 5-91, (GEN(1,J) , 1 = 1, M) 

501 FORMAT! 'O' , * SUBGENERA TOR= • , 10 1 6 ) 

PR I NT 100 1, DELTA 
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10)1 F0RMAT!'-'»'DELTA=', 13) 

PRINT6G2 

602 FORMAT! 'O' ,T2, 'BITS' ,T19,' ACCURACY' ,T37, 'NOISE COUNT: BIT1»BIT2» 

1BIT3, ETC.') 

DO 600 .1=1 , NPLOT 

600 PR[NT6-<4,KMTPLT( J) ,FPLOT(J) ,<PMPLOT( J, I) ,I=1,N) 

604 FORMAT! • • , T2 , 1 7 , T15,IH6. 8 , T39, 8 1 9 > 

PUNCH 1 105.MAXBIT 
1105 FORMAT! I7»T73, 'MAXBtT* > 

PUNCH 1 1 35 » N 

1135 FORMAT! II, T73»'N' ) 

PUNCH1 101 ? M 

11C1 FORMAT! II, T73,'M* ) 

PUNCH 1 1 36»NR0WS 

1136 FORMAT! 17, T73, 'NROWS' ) 

PUNCH1137,NCOL 

1137 FORMAT! 17, T73, 'NCOL* ) 

DO 1566 J=1,N 

1566 PUNCH1I03,J, (GEN! I,J) »I=1,M> 

1103 FORMAT (T7 3, 'GEN6R! * , 1 1, • J',Tl,9ID 
PUNCH i 107, DELTA 
1107 FORMAT! 13, T73, 'DELTA* ) 

PUNCH11 *’9, STON 

1109 FQRMAT(D15.7jT73» ' S/N ' ) 

PUNCH3000 , NXMTD 

3C0G FORMAT! 1 20, T73 , ' NXMTD ' » 

PUNCH3 lt.U, I XSQUR 

3100 FORMAT! 120, T73, ' IXSOUR' ) 

DO 1567 J=1,N 

1567 PUNCH 3200, I XNOI S ! J ) , J 

3200 FORMAT! 120, T71,* IXNOIS! ’,11,' )* ) 

PUNCH3300, ERRORS 

3300 FORMAT! I2D.T73, 'ERRORS’ ) 

PUNCH340U, NTOTAL 

3400 FORMAT! 120, T73, ’NTOTAL* ) 

DO 1569 J=1,,N 

1569 PUNCH31D1,J,!PN0ISE! IOUT+I,J) ,I=1,DELM1) 

3101 FORMAT (T73, 'PNOISE', I 1, Tit 72 II) 

PUNCH32C1, ! PRSOUR ( IOUT+I ) ,1=1, DELHI ) 

3201 FORMAT (T73, 'PRSOUR* ,T1, 7211) 

DO 1570 J= 1 » N 

1570 PUNCH350G , P ARCK ( J ) , J 

350C FORMAT! 120, T7 3, 'PARCK! • , II, ’>• ) 

DO 4000 I=1,N$TAT5 

4000 PUNCH3002, GAMMA! I ) , I 

3002 FORMAT(D23.16,T73, ' !=' , 14) 

DO 4001 1=1, NS TATE 

DO 4001 J= IDEPCP , NCOL 

4001 PUNCH 3 30 l ,I,J,SURVIV(I,J) 

3301 FORMAT! T70, * 1 = ' ,14,’ J= • , 1 2 ,T 1 , 1 20) 

PUNCH3401, » SHIFT2! I ), 1 = 1, M) 

3401 FORM AT !T73, 'SHIFT' ,T1, 7211 ) 

STOP 

END 
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TESTBT PROGRAM 


***TESTBIT*** CHECKS VALUE 

CF ANY BIT IN A FILL WCRC IN RAIN STCRAGE. 

♦ 

Rl- 

-ADCS CF ARGUMENT LIST 

* 


0111 ACDR 

OF 

HORD 

* 


Ail) AC CR 

CF 

TESTBIT 

* 


8(1) ACCR 

FOR 

RETURN CCCE 

♦ 

L 

2tCU> 


LC AD ACCR CF KCPC TC BE SHIFTED 


L 

A t A 1 1 ) 


LCAC ACCR CF TESTBIT 


l 

C «C ( A ) 


LCAC TESTBIT 


LA 

3 »C 


SET FCR CCHPARE 


L 

A f C ( 2 ) 


LCAC TEST NCRC 


LA 

5*32 


LCAC 32 FCR SUBTRACT ICN 


SR 

5tC 


FINC CCRPLEHENT 


SLL 

A t C ( 5 ) 


SHIFT SC TEST BIT IS IN THE SIGN PCSITICN 


CR 

A ♦ 3 


CCNPARE RESULTS AGAINST ZERO 


8M 

NmEG 




LA 

3*1 


LCAC CNE IF NEGATIVE 

KCTNEG 

ECU 

* 




L 

5,em 


LCAC ACCR CF RETURN WCRC 


ST 

3 » C ( 5 ) 


STORE BIT IN RETURN WCRC 


STOP 

END 
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CCPYAR PROGRAM 

***C0PYARAY***CCP1E$ AN ARRAY THROUGH THE USE CF THE MVCL INSTRUCTION. 

* INPUT: 

* R1 — ADCP OF ARGUMENT LIST 

* 0(1) ACER OF INPUT ARRAY 

* 4(1) AECR OF OUTPUT ARRAY 

* 0(1) ACER OF NUMBER OF RONS IN THE ARRAYS 

* 12(1) ACCR CF NUMBER CF COLUMNS IN THE ARRAYS 

* 

CCPYARAY ECU * 

L 2,4(1) ACCR CF ARRAYC 

L 4 , C ( 1 ) ACCR CF ARRAY I 

L 6,0(1) LCAC ACCR CF NUMBER CF RCWS 

L 5,12(1) LCAC ACCR CF NUMBER OF COLUMNS 

L 5,0(5) LCAC NUMBER CF COLUMNS 

MH 5,2(6) MULTIPLY NUMBER CF ROWS BY COLUMNS 

SLA 5,2 (RONS 4 COLUMNS) * 4 

LR 3,5 COPY LENGTH IN R3 FOR MVCL 

MVCL 2,4 COPY ARRAY 

STCP 

ENC 
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SFIfT PRCGRAK 

***CCPYRCNS***CCNTA INS THE LOGIC NECESSARY TC CCPY ANC SHIFT 

* THE SPECIFIC RCWS AS SPECIFIED. CCPYRCWS HILL CCPY RCWI 

* INTO RC WC AND ROWl, SHIFT THESE ROWS TC THE LEFT CNE eiT, 

* ANC SET THE RIGHTMOST 6 IT OF RCWO ANC ROWl TC ZERO AND CNE 

* RESPECTIVELY. 

* INPUT : 

* R 1 — ADCR OF ARGUMENT LIST 

* 0 < 1 > ACOR OF INPUT ARRAY ( ARRAY I ) 

* A<1) ACDR OF OUTPUT ARRAY I ARRAYC ) 

* 8(1) ACER OF NUMBER OF RCWS 

* 12(1) ACOR OF NUMBER CF COLUMNS 

* 16(1) ADCR CF ROM I 

* 2C( 1 ) ADCR CF ROMO 

* 24(1) ADCR OF ROWl 


CCPYRCWS EQU 
LA 
L 
L 

SLA 

L 

L 

L 

L 

LR 

BCTR 

L 

MH 

LR 

LR 

LR 

L 

L 

L 

L 

L 

L 

BCTR 

BCTR 

BCTR 

AR 

AR 

AR 

SLA 

SLA 

SLA 

AR 

AR 

AR 

MVC 

LA 

LR 

BAL 

MVC 

LA 

LR 

BAL 

BCTR 

CR 

BE 


7 » C 

12 . 8 ( 1 ) 

12 . 0 < 12 ) 

12,2 
6,0 (1) 
6 , 12 ( 1 ) 

8 , 0 (8 ) 

2,4(1) 

3.8 

3,0 

4 * 6(11 

3,2(4) 

5.3 
1C, 3 

11.3 
3,16(1) 
4,20(1) 
5,24(1) 
3,0(3) 

4, C (4) 

5, C<5> 

3.0 
4 ,C 
5,C 

9.3 
1C, 4 
11,5 

5.2 

10.2 
11,2 

5.6 
10,2 
11,2 

0(4,10,0(5) 

C,C 

2,10 

14, SLICE 
0(4, ID, 0(5) 
0,1 
2,11 

14, SHOE 

8.0 

8.7 
THRU 


ZERO CUT R7 

UO AC ADCR CF NUMBER CF ROWS 
LOAD NUMBER CF PCWS 
MULTIPLY BY FCUP FCR DISPLACEMENT 
LOAD ADCR CF ARRAY I 
LC AC ACCR CF NUMBER CF CGLUMNS 
LOAD NUMBER CF COLUMNS 
LC AC ACCR CF ARRAYC 
COPY NUMBER CF COLUMNS FROM R8 
SUBTRACT CNE FPCM NUMBER CF COLUMNS 
LCAC ACCR CF NUMBER CF ROWS 
COMPUTE RCWSMCCLUMNS-l) 

#** 

* * CCPY RCWS* (COLUMNS*!) 


* LCAC ACDRS CF RCW I »RCWC , ANC RCN1 


* LCAC RCW I , RCWO , AND RCW1 


* DECREMENT EACH BY CNE 


* ACC PRODUCT ANC ROWS MINUS CNE 


* MULTIPLY BY 4 TC OBTAIN DISPLACMENT 


♦ * COMPUTE ACTUAL ADDRESSES 

*** 

CCPY LAST WCRC CF RCWI TC RCWO 
INPUT FCR RCWO 
INPUT FCR SLICE 
BRANCH TC SLIDE 

CCPY LAST WCRC CF RCWI TC PCW1 

INPUT FCR. ROWl 

INPUT FCR SLICE 

BRANCH TC SLIDE 

DECREMENT COLUMN CCUNT BY CNE 

COMPARE COLUMN COUNT AGAINST ZERO 

IF ECUAL ENC 
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PART C 



OECACCR 

ECU 

* 



SR 

9*12 

*** 


SR 

10,12 

* * DECREMENT ACDRS 0Y RCW WEBBER 


SR 

11*12 

*** 


MVC 

0(4*10), 0(9) 

CCPY NEXT RCWI KEPBER TC RCWO 


LR 

2*10 

SET UP FCR SLICE— CARRY BIT ALREACY 


0 A L 

14, SLICE 

BRANCH TC SLICE 


PVC 

0(4,11) ,0(10) 

CCPY RCWO PEHBER INTO RCW1 PEMBEP 


8CT 

8*CEC ADCR 

DECREMENT ANC TEST CCLUPN CCUNT 

THRU 

ECU 

STCP 

♦ 


ARGLSTAD 

CS 

F 



CIST 


***SLICE+**SHIFTS 

* SETTING THE 

* AECR CF THE 

* CVER BIT IS 

* 


CNE FULL WORD IN' MAIN STORAGE TC THE LEFT ONE eiT 
RIGHTMOST BIT IN THE WCRC AS INCICATEO BY RC. THE 
WORD TO BE SHIFTED IS CCNTA INEC IN R2 AND THE CARRY 
PLACEC IN RC. 


SLICE 

ECU 

* 



STH 

3, A, SAVEREGS 

SAVE WORKING REGISTERS 


LA 

A,C 

LOAD ZERC FCR CCPPARE 


L 

3.CC2I 

COPY WCRC TC BE SHIFTED 


CR 

3 * A 

CHECK FCR CARRY 


SLL 

3,1 

SHIFT CVER CNE BIT 


8NL 

CAPRYG 

BRANCH ACCCRCINGLY 

CARRY 1 

EQU 

* 



LA 

A, 1 

WORD NEGATIVE* THEREFORE SET CARRY 

CARPYO 

AR 

3,0 

SHIFT IN BIT 0 Y ACCITICN 


LR 

0 , A 

STORE CARRY VALUE 


ST 

3,0(2) 

REPLACE SHIFTED WCRC 


LM 

3 , A , SAVEREGS 

RESTORE WORKING REGISTERS 


BR 

IA 

RETURN TC CALL— 0AL 14*SL1DE 

SAVEREGS 

CS 

2F 



END 




BIT 
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SHORT CONSTRAINT LENGTH RATE 1/2 'QUICK-LOOK 1 CODES 


l. Introduction: 

A binary convolutional code of constraint length K and rate R = £ is com- 
pletely specified by a set of two generators which in transform notation have 
the form 


0 ) 


(D) = g Q (j) + g t (j) D + g 2 (j) D 2 + . . . + 9 K _ 1 (j) D K_1 (j = 1,2) 


with coefficients from GF(2). (Throughout we assume the codes are nondegener- 
ate, i.e., at least one of 9 q^ and 9g^ are at least one of g^^^ and 

( 2 ) \ 
g ^_1 are one) . I f 

I (D) = i Q + i t D + I 2 D 2 + . . . 


is a sequence of binary information digits, then the result of applying I (D) 
to the encoder is 

T^ (D) = 1(D) (D) = + tj^ D + t 2 ^ D 2 + . . . (j = 1, 2) 


so that for each information digit i^ the encoder produces a block of two digits 

r (1) (2)n 

Lt k v , t^ J that are functions of i^ and the previous K-1 information digits. 

The linear sequential circuit that performs this operation consists of a shift 
register whose K stages are connected to two modulo-2 adders in accordance with 
the coefficients of G^ (D) and G^ (D) , respectively. The outputs of the 
adders at time k then constitute the block [t^^ , t^ 2 ^]. For convenience we 
denote the sequence of these blocks by T(D). 

In certain situations such as system check-out it is desirable to be able 
to recover the information sequence from the encoded sequence. Massey and Sain 
(1968) have shown that this is possible if and only if the code is noncatastroph ic , 
i.e-, if and only if 

gcd [G^ (D), (D)] = D^ 
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for some Z k 0. In this case, there always exists a linear sequential circuit 
that produces I (0) with a delay of exactly L digits for any integer L > £, and 
it is completely described by two generator polynomials P^ (D) and (D) 

that satisfy 

(D) G^ (D) + P^ (D) (D) = D £ 

To illustrate these ideas we consider the code 
G^ (D) = 1 + D + D 2 + D 3 + 

G^ (d) - 1 + D 2 + D 3 + D 5 + D 6 

This code has a constraint length K = 7 and its circuit realization is shown 
in Figure 51. If the input sequence is 

1 0 0 1 1 01 0 1 ... 
then T^ (D) and T^ (D) are given by 

1 1 1 0 0 0 0 0 0 ... 


and 

1 0 1 0 1 0 0 1 1 ... 

respectively, and the encoder output sequence will be 

1 1 101 1000100000101 ... 

Since G^ (D) and (d) are relatively prime, an inverse circuit with delay 

( 1 ) ( 2 ) 

zero exists and we may easily prove that P v 1 (D) and P v ' (D) are given by 

P^ (D) =1 + D + D 2 + 

P^ (D) = D 2 + 

Two versions of the circuit realization are shown in Figure 52. 

Suppose now that the encoder output sequence T (D) is transmitted over a 
noisy channel prior to its inversion. Then, of course, the resulting sequence 
I (D) will generally not be a perfect match of the original information sequence 
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I (D) . In fact, Massey and Costello (1971) have shown that over the binary sym- 
metric channel and at high signal-to-noi se ratios the probability of an error in 

A 

1(D) is related to the probability of error in the channel by 

P? = Ap BSC 

where A is the error amplification factor given by 

A = W[P^ (D)] + W[p^ (D)] 
and W[P^ (p)] denotes the Hamming weight of P^ (p). 



FIGURE 51. ENCODING CIRCUIT FOR THE CODE 
11110 0 1 
10 110 11 


In our example above A has the value 7, so that an error in 1 (d) is seven 
times more likely than an error in the channel. This is quite obvious from 
Figure 5^h. For a single error in the channel will, as it propogates through 
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T(D) 


FIGURE 52b. ALTERNATE INVERSE CIRCUIT FOR THE CODE 



11110 0 1 


Sk 


1 0 


1 


1 


0 1 1 




the circuit, produce 7 errors in the output of the adder, assuming that the 
channel errors are spaced far enough apart. 

For low signal-to-noise ratios the simple reasoning leading to (1) no 
longer applies and the value of the error amplification must be determined 
empirically. Figure 53 shows the result for the code in the above example. 

Consider next the system configuration of Figure 5^. 

At high signal-to-noise ratios a well designed decoder will be able to 
correct the overwhelming majority of the errors introduced in the channel and 
deliver an essentially perfect copy of 1(D). If we then compare this output 
with that of the encoder inverse we obtain an indication of the signal-to-noise 
ratio in the channel . 

With the binary symmetric channel, for example, we can get a good estimate 
of Pj by computing the ratio of the number of ones in which the outputs of the 
decoder and the encoder inverse differ to the total number of digits processed. 
Using (l) we are then able to determine the value of P Dcr . 

The surprising fact is that this scheme also works for low signal-to-noise 
ratios, where the decoder output also includes errors, and produces a one-to-one 


relationship between p R _ r and the measured quantity, which we denote by p_ . 


Figure 55 shows the simulation results for the code In our previous example, 


the binary symmetric channel and a 32 bit path length Viterbi decoder. 


From Figure 53 It is clear that if one attempts to reconstruct the original 
information sequence at the channel output without benefit of decoding, it is 
desirable to have a code with as low a value of error amplification as possible. 
The best in this regard are the socalled systematic codes for which one of the 
(D) is one and the other equals zero, resulting in A = 1. Unfortunately, 
the error correcting capability of these codes is markedly inferior to that of 


certain nonsystematic codes when used in conjunction with sequential or maximum 


likelihood decoding algorithms. 


S5 



ERROR AMPLIFICATION FACTOR 









BSC 




For nonsys temat i c codes the lowest possible value of A is 2 and is attained 
by the socalled quick-look codes [2], Our purpose in this note is to investi- 
gate their relevant characteristics and in the process we obtain a number of in- 
teresting and practically useful results. Since our primary motivation is the 
application of quick-look codes to Viterbi decoding, we restrict consideration 
to constraint lengths less than eight. 

II, Quick-Look Codes: 

We define a rate i quick-look code as any code in which the two generators 
differ In exactly one coefficient. Then 

G (1) (D) + G (2) (D) - D L 

for some 0 < L < K - 1 and an inverse circuit with delay L and error amplifica- 
tion factor A = 2 is given by 

P^ (D) = P’ 2 ^ (D) = 1 

This, of course, amounts to nothing more than the modulo-2 addition of 
T (1) (D) and T^(d). Hence the word 'Quick-Look 1 [2]. 

Since we are dealing with nondegenerate codes only, it follows easily that 
all quick- look codes have 

gcd [G (,) (D), G (2) ([»] - 1 

Thus, there always exists an inverse with delay zero, which is generally dif- 
ferent from the quick-look inverse if l > 0. 

For example, when L = 1 , the zero delay inverse takes the form 

p (l) (d) = ( - D) (i t j) 

and its error amplification factor at high s ignal-to-noi se ratios is 
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i ? j 


A = W[G^ (D)] + W[G^ (D)] - 2 
For L = 2 the zero delay inverse becomes 

/,) 1 + (1 + a D) G (j) (D) 

P U) (D) = 1 

D 2 

Here A has the same value as above if = 0 and is a function of the coefficients 
of G^ (D) and G^ (D) If a^ = 1 . 

As a concrete example, consider the constraint length 5 code 

(D) = 1 + D + D 2 + 

(D) = 1 + D + D^ 

Clearly, L = 2 and the quick-look inverse circuit takes either of the forms in 
Figure 56 . 

The inverse circuit with zero delay is given by 
P^ (D) = 1 + D 2 + D 3 

P^ (D) = D + D 2 + D 3 

and Figure 57 shows the two alternate configurations for this case. Mote that 
the error amplification factor increases from 2 to 6 over the quick-look 
inverse. 

ill. Maximum Free Distance Quick-Look Codes: 

One commonly accepted measure of the performance of a convolutional code 
in conjunction with sequential or maximum likelihood decoding algorithms is 
free distance. For the codes under consideration here this is simply the smallest 
nonzero number of ones in the set of semi - i nf I n i te output sequences of the 
encoder . 

Our objective is to find quick-look codes of constraint lengths 3 s K < 7, 
with as large a free distance as possible. 
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I (D) = DIGITS 5 , 7, 9, . . . 
CIRCUITS FOR THE CODE 
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110 0 1 


To narrow the search for such codes we first note that the maximum free 
distance of any rate noncatastrophtc convolutional code is bound by 

/K+2; 3sKs6 

d f < 

(K + 3; K = 7 

and that there always exists a code for which equality holds (Larsen, 1973). 

Second, since the input sequence 100 . . . produces as output sequence 
from each modulo-2 adder of the encoder the coefficients of the respective 
generator polynomial, the free distance of any code is evidently bounded by 

d f < W[G {1) (D)] + W[G (2) (D)] 

Finally, if G*(o) denotes the reciprocal polynomial of G(D), then the codes 
(D) , G^ (D) 
and 

G^*(D), G^*(D) 
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are equivalent. This follows readily from the relation 
[I (D) G (D)]* = I* (D) G* (D) 

and the fact that the weights of a polynomial and its reciprocal are identical. 

We can, therefore, restrict our search to quick-look codes with delay 
L < [K/2] and an appropriate number of ones in the generator polynomials (the 
square brackets denote the integer part). 

Now let L = 0. Then each 1 in the input sequence 1(D) will produce a 1 
in the output sequence T(D) as it enters the encoder shift register and in ad- 
dition the last 1 in l(D) will produce two I's in T (D) as it enters the last 
stage of the encoder. Therefore, 

W[T CD) ] 2 2 + W[l (D)] 

and it follows that in testing whether a code has free distance less than 
only input sequences with fewer than d^ - 2 ones need to be considered. 

Since Bah] and Jelinek (1971) have shown that without loss of generality 
input sequences with zero-runs of length K - 2 or more may likewise be ignored, 
it follows that the length of the Input sequences that must be tested does not 
exceed 

(d f -4) (K - 2) + 1 

For L > 0, the first 1 in 1(D) produces two I's in T (D) as It enters the 

encoder and another 1 as it enters the (L +.l)st stage of the encoder. Every 

subsequent 1 in 1(D) likewise produces a 1 in T(D) as it enters the (L + l)st 
stage. In addition, the last 1 in T (d) results in two I's in T (D) as it 

enters the last stage of the encoder. Thus, the total number of ones in the 

output sequence satisfies 

W[T (D)] > 4 + W[l (D)] 
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and we can restrict consideration to input sequences with fewer than - 4 
ones and length no larger than 
(d f - 6) (K - 2) + 1 

Using these principles we tested all quick-look codes of constraint length 
3 s K s 7. Table 3 summarizes our results. Note that for 3 < K < 6 the 
best quick-look codes are comparable to the best genera] nonsystematic codes, 
whereas for K = 7 the free distance of the best quick-look codes is one less 
than the maximum achievable. 

We also remark that the quick-look codes with L = 0 are uniformly inferior 
to those with L > 0, a result that reinforces the notion that among the best 
codes of a class there is always one whose generators possess complementarity 
(Bah 1 and Jelinek, 1972). 

Although under normal circumstances free distance is a good indicator of 
a code's error correcting capability, this measure nevertheless depends only 
on the code and thus completely ignores the nature of the channel and the 
decoding algorithm. Even with the channel and decoder fixed, differences in 
the weight spectra of two codes with the same free distance can give rise to 
different decoder bit error rates. 

For these reasons we have computed the decoder bit error rates of selected 
codes from Table 3 used over the binary symmetric channel and in conjunction 
with a Viterbi maximum likelihood decoding algorithm of 32 bit decoder path 
lengths. The results are presented in Figure 58 . Note that these quick-look 
codes compare favorably to the best nonsystematic codes obtained in [7] and 
the complementary codes given by Jelinek and Bahl (I 969 ). 

in Figure 59 we show the error amplification factor for the same set of 
codes as above, as a function of the s i gna 1 -to-noi se ratio of a binary symmetric 
channel . 
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K, 

Code § 

(octal) 

L 

d f 

d fmax 

3 

1 

7 

1 

5 

5 

4 

2 

17 

1 

6 

6 

5 

3 

33 

1 

7 

7 

4 

35 

2 

7 

7 

6 

5 

67 

1 

8 

8 

6 

75 

1 

8 

8 

7 

7 

153 

1 

9 

10 

8 

163 

1 

9 

10 

9 

127 

2 

9 

10 

10 

135 

2 

9 

10 

11 

165 

2 

9 

10 


12 

171 

2 

9 

10 

mm 

175 

2 

9 

10 

14 

133 

3 

9 

10 


Best Rate 1/2 Quick-Look Codes 
TABLE 3 


Finally, Figure 60 presents the relationship between actual and measured 
channel bit error rates for the same codes, the binary symmetric channel and a 
32 bit path length Viterbi decoder. 
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FIGURE 58. VITERBI MAXIMUM LIKELIHOOD DECODING 
PERFORMANCE A = 32 
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FIGURE 59. ERROR AMPLIFICATION FACTOR FOR CODES #1, 2, 3, 7 
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FIGURE 60. MEASURED VERSUS ACTUAL CHANNEL BIT ERROR 
PROBABILITY FOR CODES #1, 2, 3, 7 
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